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Introduction 


Local  area  networks  (LANs)  are  designed  to  support  high  bandwidth  communica¬ 
tions  among  a  large  number  of  users  within  a  local  geographical  area.  When  a  single 
channel  (coaxial  cable,  fiber  optic,  or  radio  multi-access  channel)  is  available, 
and  the  users  are  allowed  to  access  the  channel  randomly  (whenever  they  have  some¬ 
thing  to  transmit),  a  distributed  control  random  access  algorithm  (RAA)  may  be  used 
to  coordinate  the  overall  user  activity.  Wien  the  end-to-end  propagation  delay  of 
the  LAN  is  small,  as  compared  to  the  transmission  time  of  a  packet,  then  the  users 
can  determine  the  channel  activity  in  a  short  amount  of  time,  through  "channel 
sensing"  operations. 

The  most  well  known  RAAs  for  LANs  belong  to  the  class  of  "Aloha-type"  algorithms, 
like  the  Non-persistent  CSMA,  CSMA/CD  and  their  variations  [1-4].  Although  the  above 
algorithms  are  easy  to  implement,  they  have  inherent  long-term  stability  problems. 

The  recognition  of  this  Instability,  by  some  researchers,  led  to  the  "Aloha-type 
with  retrnasmission  control"  [5]  and  the  "Tree-Search-type"  [6-12]  classes  of 
algorithms.  The  algorithms  in  those  two  classes  are  stable,  but  they  are  inherently 
characterized  by  a  major  operational  drawback.  An  integral  part  in  their  operation 
is  the  assumption  that  each  user  senses  the  channel  continuously,  even  when  he  is 
not  yet  part  of  the  system.  This  full  channel  sensing  assumption  equivalently  implies 
that  each  user  "knows"  precisely,  and  at  all  times,  the  lag  between  the  time  instant, 
when  the  system  started  operating,  and  the  time  instant  of  his  own  arrival  in  the 
system.  This  extreme  synchronization  is  clearly  nonfeasible  in  most  real  systems, 
where  activation  of  new  users,  user  mobility,  and  actions  of  higher  level  protocols 
discarding  packets  already  in  the  system  [10,11],  disrupt  the  feedback  sensing 
continuity.  In  summary,  the  "Aloha-type"  algorithms  have  poor  performance  due  to 
Instabilities,  and  the  "Aloha-type  with  retransmission  control"  and  the  "Tree- 
Search-type"  algorithms  are  nonimplementable,  in  most  systems,  due  to  extreme 
synchronization  requirements. 


In  this  paper,  we  consider  the  class  of  "Limited  Channel  Sensing"  algorithms, 
[13-16,  21].  The  algorithms  in  this  class  are  stable,  and  they  can  attain  high 
quality  performance.  In  contrast  to  the  "Tree-Search-type"  algorithms,  the  "limited 
channel  sensing"  algorithms  require  that  each  user  sense  the  channel,  only  from  the 
time  that  he  generates  a  packet,  to  the  time  that  this  packet  is  successfully  trans¬ 
mitted  (limited  channel  sensing) .  This  level  of  channel  sensing  is  clearly  feasible 
in  most  LANs.  We  propose  and  analyze  "Limited  Channel  Sensing"  algorithms  for  the 
CSMA  and  CSMA-CD  channels.  In  contrast  to  the  algorithm  in  [18],  the  proposed 
algorithms  operate  asynchronously.  That  is,  we  assume  that  the  users  can  not  dis¬ 
tinguish  slot  boundaries;  thus,  packet  transmissions  start  at  random  times,  rather 
than  at  the  beginning  of  some  slot.  Asynchronous  algorithms  are  of  greater 
practical  importance  than  synchronous  algorithms,  especially  for  local  area  networks 
since  there  is  no  need  for  the  users  to  maintain  a  global  time  base.  In  addition, 
such  algorithms  may  out  perform  their  synchronous  counterparts,  in  some  important 
local  area  environments  (as  we  will  see  later) . 

Utilizing  the  regenerative  character  of  the  stochastic  processes  that  are 
associated  with  the  random  access  channel,  we  derive- lower  bounds  on  the  maximum 
stable  throughput,  and  tight  upper  and  lower  bounds  on  the  induced  mean  packet 
delay.  These  results  indicate  that  the  proposed  algorithms  are  inherently  stable, 
and  that  they  combine  good  performance  with  limited  channel  sensing  and  modest 
synchronization  requirements. 

2 .  User  and  Channel  Model 

We  assume  that  an  infinite  population  of  identical,  independent,  bursty,  packet 
transmitting  users  share  a  common  communication  channel.  We  model  the  cumulative 
packet  arrival  process  as  homogeneous  Poisson,  with  intensity  \  packets  per  unit  of 
time.  We  conveniently  take  the  packet  length  to  correspond  to  our  unit  of  time.  We 
also  assume  that  the  propagation  delay  between  any  two  users  in  our  system  is  at 
most  a,  where  a«l.  We  note  that  the  above  Poisson  user  model  best  represents 


environments,  such  as  the  packet  radio  and  the  Ethernet,  where  bursty  independent 
users  enter  and  depart  the  system  randomly. 

We  consider  limited  channel  sensing  and  ternary  feedback.  That  is,  each  user 
senses  the  channel  continuously,  from  the  time  Instant  when  he  generates  a  packet, 
to  the  time  instant  when  this  packet  is  successfully  transmitted,  and  he  can  dis¬ 
tinguish  without  error  among  the  following  channel  states:  a)  idle  (no  transmission) 
b)  success  (transmission  of  a  single  packet)  c)  collision  (simultaneous  transmission 
of  at  least  two  packets) .  We  assume  that  a  collision  results  in  complete  loss  of 
the  information  included  in  all  the  involved  packets;  thus,  retransmission  is  then 
necessary.  We  assume  that  the  users  can  not  determine  slot  boundaries.  We  thus 
consider  unslotted  channel,  and  asynchronous  transmissions.  We  initially  introduce 
the  following  assumptions: 

Al)  The  propagation  delay  between  any  two  users  in  the  network  is  exactly  a. 

A 2)  A  collision  is  distinguished  from  a  successful  transmission,  by  all  the 
users  in  the  system  (those  who  are  then  sensing  the  channel),  if  it  lasts 
Y  time  units. 

We  point  out  that  assumption  Al  is  initially  adopted  for  analytical  reasons. 

In  practice,  the  propagation  delay  between  any  two  users  is  a  random  variable, 
depending  on  the  relative  position  of  those  two  users  in  the  network.  We  will 
discuss  the  relaxation  of  assumption  Al  in  section  6.  Assumption  A2  is  realistic, 
where  the  length  y  is  a  network  characteristic  determined  by  the  encoding  of  the 
packets  in  baseband  transmission,  or  by  the  modulation  technique  used  in  broadband 
transmission.  Note  that  assumption  A2  implies  that,  in  addition  to  transmitting 
users,  non-transmitting  users  also  have  the  capability  to  detect  collisions,  provided 
that  they  are  sensing  the  channel  during  the  collision;  (receive  mode  collision 
detection,  IEEE  Standards  Committee,  Project  802,  [25]). 

3.  An  Algorithm 

We  consider  the  user  and  channel  model  in  section  2,  and  we  adopt  assumptions 
Al  and  A2.  We  assume  that  all  users  know  then  the  common  propagation  delay,  a. 


and  the  common  parameter,  y,  and  we  describe  an  asynchronous  limited  channel  sensing 
algorithm  implemented  by  each  user  in  a  distributed  fashion.  We  name  the  algorithm, 
SALCS. 

Each  user  has  in  memory  the  parameters,  a  and  y,  and  two  common  parameters,  m 
and  n,  that  are  both  positive  integers,  subject  to  optimization  for  throughput 
maximization.  The  time  instant  when  the  user  generates  a  new  packet,  he  immediately 
starts  sensing  the  channel,  and  he  simultaneously  initializes  the  SALCS.  He  maintains 
continuous  channel  sensing  until  the  successful  transmission  of  his  packet;  upon  the 
occurrence  of  this  event,  he  stops  sensing  the  channel.  For  the  implementation  of 
the  SALCS,  the  user  has  a  counter,  whose  value  of  time  t  is  denoted  r^.  The  values 
of  the  counter  dictate  the  operation  of  the  algorithm,  which  is  described  as  follows: 

A.  Transmission  rule 

The  user  transmits  whenever  his  counter  value  equals  "0". 

B.  Initial  counter  value 

Let  the  user  generate  a  new  packet  at  time  to»  Then,  he  immediately  starts 
sensing  the  channel,  and  he  sets: 

!0  ;  if  he  senses  the  channel  idle  at  t  . 

° 

-1  ;  if  he  senses  the  channel  busy  at  tQ. 

C.  Special  action  after  collision 

Let  at  time  t  the  user  start  a  transmission.  Let  t'  be  the  first  time  instant 
after  t  that  the  user  senses  the  channel  state  as  collision,  while  still  transmitting. 
Then, 

a)  If  t'-t  >_  y  +  a,  the  user  extends  his  transmission  to  the  time  instant 
t"=t+y+2a.  At  t",  he  aborts  his  packet,  and  he  considers  t"  as  the  time 
instant  when  the  collision  ends. 

b)  If  t'-t  <  y  +  a,  the  user  aborts  his  packet  at  t',  and  considers  t'  as 
the  end  of  the  collision. 

D.  Recurs ions  of  counter  value 


Let  at  time  t  the  user  sense  the  channel  busy,  and  let  his  counter  value  be 


then 

idle 


rt<  Let  t'  be  the  first  time  instant  after  t,  that  the  user  senses  the  channel 
Then, 

Rule  1:  The  counter  value,  r^,  remains  unchanged,  until  the  time  instant  t* . 
Rule  2:  At  time  t',  the  counter  value  is  updated  as  follows. 

a)  If  rt  =  -1,  then, 

rt<  =  M  ;  0  <_  M  £  m-1,  with  probability  m  ^ . 

b)  If  at  (t1)  the  user  senses  the  channel  state  as  success,  and 

r  >_  1,  then. 


r 


t' 


r  +  m-1 


If  r  =  0,  the  user  has  been  successfully  transmitted  at  t'. 


c)  If  at  (t')  the  user  senses  the  channel  state  as  collision,  and 
rt  0,  then. 


V 


rt+m+-n-l 


if  r  >  1 


-1 


|  m-l+N;lj£lKn,w.p.  n  ;  if  r£  =  0,  and  the  user  acted 

as  in  C.  b)  upon  sensing  the 
collision. 


-1 


if  r  =  0,  and  the  user  acted 

as  in  C.  a)  upon  sensing  the 
collision. 


Rule  3:  After  time  t',  the  counter  value  is  updated  as  follows. 


a)  If  r£i  >  1,  the  user  decrements  his  counter  value  by  one,  every 
a  time  units,  after  t',  to  a  value  not  less  than  0,  until  the 
time  instant  t",  that  he  first  senses  the  channel  busy  again. 
Then,  he  maintains  the  counter  value  rt«  unchanged,  until  he 
first  senses  the  channel  idle  again.  He  then  repeats  rule  2. 

b)  If  r£  =  0  and  rt,  =  -1,  then. 


r  ■■  -1  ;  for  t'  <  T  <  t'  +  a 


r 


“  m-l+N, 


1  <  N  <  n. 


with  probability  n 


-1 


f- f  4-rt 
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The  user  decrements  his  counter  value  by  one,  every  a  time 
units,  after  t'-fa,  to  a  value  not  less  than  0,  until  the 
time  instant  t" ,  that  he  first  senses  the  channel  busy  again. 
Then,  he  maintains  the  value  rt"  unchanged,  until  he  first 
senses  the  channel  idle  again.  He  then  repeats  rule  2. 

We  note  that  in  the  SALCS,  packets  do  not  use  the  distinction  between  success 
and  collision,  until  they  first  sense  an  idle  channel  state.  This  property  actually 
represents  an  advantage,  since  this  distinction  may  be  impossible  when  the  first 
sensed  idle  channel  state  appears  shortly  after  arrival.  From  the  operation  of  the 
algorithm,  we  conclude  that  the  interval,  [t't'+a),  in  part  b)  of  rule  3  is 
necessarily  idle.  Also,  due  to  assumption  A2  in  section  2,  in  conjuction  with  the 
SALCS  operation,  we  conclude  that  a  successful  transmission  lasts  one  time  unit, 
while  a  collision  lasts,  6  =  y  +  2a,  time  units.  The  relationship  between  8,  Y,  and 
a  is  exhibited  in  figure  1. 

The  general  operation  of  the  SALCS  is  perhaps  better  illustrated,  by  introducing 
the  concept  of  a  "stack".  A  stack  is  an  abstract  storage  device,  consisting  of  an 
infinite  number  of  cells,  labelled  -1,  0,1, 2, 3,...  (figure  2).  The  number  of  packets 
that  a  cell  can  accommodate  is  unrestricted.  At  each  time  t  during  the  operation  of 
the  algorithm,  users  with  counter  value  r  can  be  thought  of  as  having  stored  their 
packets  in  cell  "r"  of  the  stack.  A  packet  is  transmitted  whenever  it  enters  cell 
"0"  of  the  stack.  Packets  are  eventually  successfully  transmitted  after  moving 
through  the  cells  of  the  stack  in  accordance  with  the  algorithmic  rules  described 
above . 

4 .  SALCS  Renewal  Properties.  Stability 

Given  the  user  and  channel  model  in  section  2,  and  adopting  assumptions  A1  and 
A2,  the  SALCS  operates  in  sessions.  The  first  session  begins  at  time  zero,  when 
the  system  starts  operating.  The  ith  session  begins  at  time  R^,  and  ends  at  time 
^i+1’  w^en  t*ie  (i+Dth  session  begins;  thus.  R^  =  0.  To  define  the  renewal  time 
instants  R j ,  J  1,  we  use  the  concept  of  a  marker,  that  operates  on  the  stack  in 
figure  2.  The  marker  is  perceived  as  an  outside  observer,  who  knows  the  rules  of 
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the  SALCS,  and  who  senses  the  channel  continuously,  with  delay  a.  At  time  R^  =  0, 
the  marker  is  placed  at  cell  til.  The  marker  stays  at  cell  til  until  time  t' ,  when 
the  channel  changes  state  from  busy  to  idle  for  the  first  time  after  R^.  From  then 
on  the  marker  moves  through  the  cells  of  the  stack  according  to  the  rules  of  the 
algorithm.  At  some  time  instant  t  t'  (before  the  session  ends),  the  marker  is 
at  cell  ti  k  of  the  stack,  if  from  the  algorithmic  operations,  in  conjuction  with  the 
marker's  channel  sensing,  it  is  deduced  that  cells  #j,  j  >_  k  are  necessarily  empty 
at  time  t.  R ^  occurs  at  some  time  instant  t",  if  at  time  t"+a  the  marker  is  placed 
at  cell  tiO.  At  R ^  the  marker  is  automatically  reset  at  cell  #1,  and  the  above 
process  is  repeated  to  determine  R^.  This  process  continues  indefinitely.  Let  us 
denote  by  L^,  the  length  of  the  ith  session.  From  the  above  description  of  the 
sessions,  in  conjuction  with  the  memory less  property  of  the  Poisson  arrival  process, 
we  then  easily  conclude  that  {L^}  is  a  set  of  i.i.d.  random  variables.  The  sequence 
forms  then  a  renewal  process,  where. 


Rx  =  0 


n+1 


R  +  L 
n  n 


n=l , 2 , . , 


(1) 


Let  *(t)  be  a  time  function,  such  that  f(t)  =  1,  during  successful  transmission 
periods,  and  f(t)  =  0,  otherwise.  Let  denote  the  number  of  successful  trans¬ 
missions  in  the  nth  session,  and  let  us  define. 


R(t) 


f(t)dt 


o(t)  & 


S  =  E{S  }  ,  L  =  E{L  } 
n  n 


(2) 

(3) 


Given  n,  the  random  variable  Sn  depends  on  the  random  variable  L^,  but  the  pairs 

(L  ,  S  ),  n>l  are  i.i.d.  Also,  due  to  (2)  we  have  S  ■  R(R  )  -  R(R  ,),  and  we 
n  n  —  n  n  n-1 

can  express  the  following  theorem  (see  [17,  thr.  3.6.1]). 


Theorem  1  If  S  <  00  and  L  <  00 ,  then  there  exists  a  real  number  o,  such  that 

£im  a(t)  =  £im  ^  -  =  y-  =  a  with  probability  1;  the  quantity  0  is  the  channel's 

t-HJO  £-*»  c  L 

output  rate. 

Directly  from  the  definition  of  a  session,  we  conclude  that  all  the  arrivals 
generated  while  the  session  is  in  progress  are  successfully  transmitted  by  its  end. 
Thus,  also  represents  the  number  of  arrivals  during  the  nth  session.  Let  X^ 
denote  the  first,  after  the  beginning  of  the  nth  session,  time  that  an  arrival 
occurs.  Let  X^,  i  >_  2,  represent  the  interarrival  time  between  the  (i-l)th  and 
the  ith  arrival,  during  the  nth  session.  Let  X  denote  the  first  time  that  an 

arrival  occurs,  after  the  end  of  the  nth  session.  Then  (see  figure  3), 

S  +1 
n 

L„ +  x  ■  L  xi  <4> 

i=l 

The  random  variables,  {X^}^^,  are  i.i.d.,  with  common  mean  A  Also,  if 
L  <  00 ,  then  S  <  °°,  where  L  and  S  are  given  by  (3).  Thus,  assuming  L  <  m,  taking 
expected  values  in  (4),  and  applying  Wald's  equation  [17],  we  obtain,  L+A  *=A(S+1), 
or. 


S  *  A  L  ;  if  L  <  “  (5) 

In  view  of  (5),  theorem  1  yields. 

Corollary  1:  If  L  <  00 ,  then  0  =  A;  that  is  the  ALCS  is  then  stable . 

Thus,  the  algorithm  maintains  the  rate,  and  it  is  thus  stable,  if  the  mean 
session  length  is  finite.  In  the  sequel,  we  investigate  the  conditions  under  which 
L  <  °°,  and  we  establish  bounds  on  L. 

a.  Mean  Session  Length.  Throughput 

For  evaluating  the  expected  value  of  a  session,  we  need  the  concept  of  a 
subsession.  Unlike  sessions,  subsessions  may  be  overlapping.  A  subsession  of 
multiplicity  k  1  starts  at  time  t  if  i)  cell  It 0  of  the  stack  is  empty  at  t  , 
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ii)  cell  It 0  of  the  stack  is  nonempty  at  time  t,  and  iii)  cell  It 0  of  the  stack 
contains  k  packets  at  time  (t4a)  .  The  end  of  a  subsession  of  multiplicity  k  is 
determined  via  a  marker.  The  marker  (as  in  the  case  of  a  session)  is  perceived  as 
an  outside  observer,  who  knows  the  rules  of  the  SACLS,  and  who  senses  the  channel 
continuously  with  delay  a.  At  the  beginning  of  a  subsession  of  multiplicity  k  ^  1, 
the  marker  is  placed  at  cell  It  1,  and  stays  there  until  the  channel  changes  state 
from  busy  to  idle  for  the  first  time  after  the  beginning  of  the  subsession.  From 
then  on  the  marker  moves  through  the  cells  of  the  stack  according  to  the  rules  of 
the  algorithm.  The  subsession  ends  at  some  time  t"  >  t  if  the  marker  is  placed  at 
cell  ItO  at  time  t"+a. 

Let  (k  1)  denote  the  expected  length  of  a  subsession  of  multiplicity  k,  and 
let  L  be  the  expected  length  of  a  session  in  (3) . 

Consider  now  figure  4.  Suppose  that  a  session  begins  at  time  0  and  that  the 
first  arrival  occurs  at  time  t.  It  is  then  easy  to  see  that. 


-Act  (Aa)k_1 
6  (k-1) ! 


(6) 


From  expression  (6)  it  is  clear  that  to  evaluate  the  expected  length,  L,  of  a 
session,  we  need  to  evaluate  the  expected  lengths  of  a  subsession  of  multiplicity 
k,  Vk>  1.  The  dynamics  of  the  SALCS  yield  the  following  proposition,  whose  proof  is 
in  the  appendix. 


Proposition  1  The  expected  lengths,  of  subsessions  satisfy  an  infinite 

dimensionality  linear  system  of  the  form. 


xk  ■  £  Vu  \ +  £k 

y=i 


k  >  1 


(7) 


;  where  the  coefficients  {A.  }and  f  are  determined  by  the  system  characteristics, 

r,  y  r 

and  are  given  in  the  appendix. 
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The  following  theorem  specifies  a  sufficient  condition  under  which  system  (7) 
has  a  solution,  with  0  £  y^  <  °°,  for  all  1  £  k  <  00 ,  which  coincides  with 

the  sequence,  of  the  expected  subsession  lengths  induced  by  the  algorithm. 

Theorem  2 

(i)  Given  a,  8,  m  £  1,  n  £  2,  there  exist  bounded  functions  of  X,  b,  b' ,  c,  c'  and 
a  positive  number  Ao(a,8,m,n)  <  1  such  that  the  system  (7)  has  a  solution,  {y^}^^, 
such  that. 


0  <  b'k-c*  £  y^  £  bk-c  ,  k  £  1 

if  X  <  X  (ot,8,m,n) . 
o 

(ii)  Given  a, 8,  given  m,n,  then,  for  every  Ae(0,Ao(a,8,m,n)) ,  =  y^,  for  all  k  £  1. 

The  coefficients  b,  b'  and  c,  c'  and  the  number  Ao(a,8,m,n)  are  derived  in  the 
proof  of  the  theorem,  which  can  be  found  in  the  appendix.  Given  the  network  parameters, 
a  and  8,  let  us  now  define, 

X.(a,8)  =  sup  Ao(a,8,m,n)  -  Ao(a,8»m*,n*)  (8) 

m,n 

Directly  from  theorem  2  in  conjuction  with  expressions  (6)  and  (8) ,  the  following 
corollary  evolves. 

Corollary  2  Given  a,  8,  the  SACLS  with  operation  parameters  m*,  n*  is  stable,  for 
Ae(0,X(a,8)) 

We  used  numerical  search  techniques  to  determine  the  optimum  values  m*  and  n*, 
which  maximize  XQ(a, 8,m,n) .  In  table  1  we  give  the  values  of  m*,  n*  and  Ma,8), 
for  representative  values  of  a  and  for  8  =  2a  and  8=1.  Given  a  and  8,  the  number 
X(a,8)  is  a  tight  lower  bound  on  the  throughput  of  the  SALCS,  subject  to  assumptions 
A1  and  A2.  Throughput  is  here  defined  as  the  maximum  Poisson  rate  maintained  by  the 
algorithm. 


5. 


Delay  Analysis 

Let  the  arriving  packets  be  labelled  n=l,2,...,  according  to  the  order  of 
their  arrival  instant.  Let  V  denote  the  delay  experienced  by  the  nth  packet  (time 
interval  between  the  time  of  its  arrival  and  the  time  of  its  successful  reception 
at  its  destination).  Also,  let  denote  the  total  number  of  packets  that  were 
successfully  transmitted  during  the  first  i  sessions.  The  sequence  *-s  a 

renewal  sequence,  since  Mq  «  0,  and  =  M^,  +  Si+1>  i  ^  0;  where,  S^,  i  ^  1, 

represents  the  number  of  packets  that  were  successfully  transmitted  during  the  ith 
session.  From  the  renewal  properties  of  the  algorihtm,  it  can  be  seen  that,  for 
every  i  ^  0,  the  process  +n)n>i»  *s  a  Probabilistic  replica  of  the  process 
{Pn)n>^.  Thus,  the  process  *s  regenerative  with  respect  to  the  renewal 

sequence  with  common  regenerative  cycle,  S,  the  number  of  packets  success¬ 

fully  transmitted  over  a  session.  The  next  theorem  is  a  combination  of  theorem  2 
and  corollary  2  of  [20]. 


Theorem  3  If  (A.l)  S  is  not  periodic,  with  S  «  E(S)  <  ®,  and  if  (A. 2) 

.  S  1  n 

T  **  E[  £  V  ]  <  00 ,  then  there  exists  a  parameter  D  such  that  D  =  £im  —  £  V  .  = 

i=l  '  ir*°°  n  i=l 

1  n 

»  tim  —  E(  £  V^) ,  with  probability  1.  Furthermore,  converges  in  distribution 

to  a  random  variable  P  .  and 

00 


D  “  E(P00)  -  T/S  <  «  (9) 


Consider  again  figure  4.  Assume  that  a  session  starts  at  time  0,  and  that  the 
first  arrival  occurs  at  time  t.  Let  denote  the  expected  cumulative  delay 
experienced  by  all  packets  associated  with  a  subsession  of  multiplicity  one,  1-Kx 
units  of  time  after  the  beginning  of  the  subsession.  Let  also  W^(k>l)  denote  the 
expected  cumulative  delay  experienced  by  all  packets  associated  with  a  subsession 
of  multiplicity  k,  $4q  units  of  time  after  the  beginning  of  the  subsession.  We 
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note  that  the  packets  associated  with  a  session  or  a  subsession  are  the  correspondingly 
successfully  transmitted  packets.  Let  us  now  define. 


T1  =  l  +  a  +  -j  +  w1 


Tk  =  0  +  a  +(k-l)  (0+  |)  +  ;  k  >  ! 


(10) 

(ID 


Clearly, 


*-E 

k=l 


-\a  ( la ) 


k-1 


(k-1) ! 


(12) 


Due  to  the  dynamics  of  the  SALCS,  and  as  in  proposition  1,  it  can  be  shown  that 
the  set,  (W  )k>^,  satisfies  a  linear  system  as  that  in  (7).  Then,  a  theorem  parallel 
to  theorem  2  can  be  expressed.  Its  proof  is  parallel  to  that  of  theorem  2,  and  it 
is  omitted. 

Theorem  4  Given  a.0,  given  m*  and  n*  as  in  (8).  let  the  number  X  (a.8.m*.n*)  ** 

~ ' — — — —  o 

**  X.(a,X)  be  as  in  (8).  Then,  for  every  X  in  (0,X(a,8))  there  exist  positive  bounded 
functions  of  X,  V^,  V^,  V^,  V^,  and  (whose  expressions  can  be  found  in  [19]), 

such  that, 

V|k2  +  V^k  +  <  Wk  <  Vxk2  +  V2k  +  V3  ;  k  >  1  (13) 


We  note  that  unlike  the  bounds  on  the  L.  's  in  theorem  2,  the  bounds  on  the  T.  's 

k  k 

are  quadratic.  Those  bounds  in  conjuction  with  (9),  (10),  (11),  (12),  (5),  (6),  and 
the  bounds  in  theorem  2,  provide  the  means  for  the  computation  of  upper  and  lower 
bounds  on  the  expected  per  packet  delay,  D.  Tighter  bounds  can  be  computed, 
however,  via  the  method  of  truncated  infinite  linear  systems,  as  in  [21].  The 
specifics  of  the  method  for  the  present  problem  is  omitted  here,  due  to  lack  of 
space,  and  it  can  be  found  in  [19].  We  only  include  the  computed,  for  various  a 
and  0  choices,  bounds.  Those  bounds  can  be  found  in  table  2,  wb»re , 
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°£  1  D  i  Du  (14) 

We  note  that  the  bounds  and  remain  extremely  tight  (actually  they  coincide 
to  the  fifth  decimal  point),  even  for  X  values  close  to  X(ct,8).  The  expected  per 
packet  delay  consists  of  two  parts.  The  expected  per  packet  access  delay,  which 
corresponds  to  the  mean  time  required  for  a  packet  to  access  the  nonnegatively 
indexed  cells  of  the  stack  (cells  It 0,1,2  ...),  and  the  expected  per  packet  contention 
delay,  which  corresponds  to  the  mean  time  elapsed  from  the  moment  a  packet  enters  the 
nonnegatively  indexed  cells  of  the  stack  until  it  is  successfully  transmitted.  The 
mean  packet  delay  characteristics  are  perhaps  better  illustrated  in  figures  5  and 
6,  where  we  plot  D £,  D  versus  X,  along  with  the  mean  packet  delays  induced  by  the 
synchronous  limited  channel  sensing  algorithm  developed  in  [18],  for  a  ■  .1,  and 
8=1  (i.e.  y  =  .8  for  a  -  .1)  or  8  ■  2a  (i.e.  Y  ■  0) .  We  point  out  that  Y  *  0 
corresponds  to  8  =  2a  for  the  SALCS,  and  it  corresponds  to  8  ■  a  for  the  algorithm 
in  [18],  while  Y  *  *8  corresponds  to  8  “  1  for  the  SALCS  and  to  0  ■  .9  for  the 
algorithm  in  [18] .  We  note  that  for  the  same  values  of  a  and  y  and  for  small  X 
values  (X  <  .25),  the  expected  per  packet  delay  Induced  by  the  SALCS  is  smaller  than 
the  expected  per  packet  delay  Induced  by  the  synchronous  algorithm  in  [18].  This  is 
due  to  the  fact  that  for  small  X  values,  the  expected  per  packet  access  delay  induced 
by  the  former  algorithm  is  smaller  than  the  same  delay  induced  by  the  latter  algorithm, 
while  the  expected  per  packet  contention  delay  is  then  almost  the  same  for  both 
algorithms.  As  X  increases,  and  the  expected  per  packet  contention  delay  becomes 
dominant  (as  compared  to  the  expected  per  packet  access  delay) ,  the  synchronous 
algorithm  out  performs  the  SACLS.  In  environments  where  assumption  A1  is  true,  but 
synchronization  is  difficult  to  achieve,  this  difference  in  performance  may  be  more 
than  off  set  by  the  ease  in  implementation  of  the  asynchronous  SALCS. 

6.  The  Relaxation  of  Assumption  A1 

In  sections  4  and  5,  we  analyzed  the  performance  of  the  SALCS,  subject  to 
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assumption  Al,  in  section  2.  That  is,  we  assumed  that  the  propagation  delay  between 
any  two  users  in  the  network  is  constant  and  identical,  and  it  equals  a  time  units. 

This  assumption  is  in  practice  true,  only  for  star  networks.  In  mobile  radio  net¬ 
works,  and  in  LANs  using  a  bus,  as  Ethernet  does,  the  propagation  delay  between  two 
arbitrary  users  is  a  random  variable.  The  value  of  this  variable  depends  on  the 
relative  distance  between  those  two  users  in  the  network,  and  it  takes  values  in 
(0,  a],  where  a  is  here  the  network  end-to-end  propagation  delay  (i.e.,  the  maximum 
propagation  delay) .  The  maximum  propagation  delay,  a,  is  a  network  characteristic, 
and  it  is  reasonable  to  assume  that  it  is  known  by  all  users.  Considering  the  SALCS 
whose  operational  parameter  is  this  maximum  propagation  delay,  a,  the  issue  is:  what 
are  its  performance  characteristics,  in  the  mobile  radio  and  the  Ethernet  environments. 

In  the  mobile  radio  environment,  let  some  user.  A,  generate  a  new  packet,  and 
start  implementing  the  SALCS  until  his  packet  is  successfully  transmitted.  The 
relative  to  A  positions  of  other  users  that  interfere  with  A's  transmissions, 
during  the  above  period,  change  then  randomly.  This  random  change  is  best  modelled 
by  random  propagation  delays.  Specifically,  the  actual  propagation  delay  between 
user  A-and  some  other  transmitting  user  is  changing  then  every  T  time  units,  and 
each  time  it  is  a  random  variable,  uniformly  distributed  in  [0,  a],  where  a  is 
the  network  end-to-end  propagation  delay.  The  period  T  is  determined  by  the  mobility 
characteristics  of  the  users;  it  is  relatively  small  for  highly  mobile  users,  and  it 
is  relatively  large  for  slowly  moving  users. 

In  the  Ethernet  environment,  each  user  selects  upon  arrival  his  position  on  the 
line,  randomly  and  uniformly.  He  remains  in  the  same  position  until  he  departs  the 
system.  If  a  is  the  line  end-to-end  propagation  delay,  it  is  equivalent  to  say 
that  the  user  selects  upon  arrival  his  propagation  delay  from  a  reference  point  (one 
end  of  the  line),  randomly  and  uniformly  in  [0,  a],  and  that  he  maintains  this  property 
until  he  departs  the  system.  Thus,  the  relative  propagation  delay  between  two  users 
in  the  system  equals  then  the  difference  between  their  propagation  delays  from  the 
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reference  point,  it  is  strictly  determined  by  the  position  of  the  users  on  the 
line,  and  it  remains  unchanged  until  the  users  depart  the  system. 

We  simulated  the  SALCS  in  the  presence  of  the  mobile  radio  and  the  Ethernet 
environments,  for  the  Poisson  user  model,  where  each  new  packet  is  a  different  user. 

We  selected  representative  values  of  the  parameter  y,  and  of  the  maximum  end-to-end 
network  propagation  delay;  we  denoted  the  latter,  a.  In  both  the  mobile  radio  and 
the  Ethernet  simulations,  we  selected  the  maximum  propagation  delay,  a,  and  a/2,  as  the 
SALCS  operational  parameter.  We  assumed  that  the  position  of  each  user  remains  practically 
unchanged,  from  the  time  he  generates  a  packet  to  the  time  this  packet  is  success¬ 
fully  transmitted.  This  latter  assumption  reflects  mobile  radio  LANs  accurately, 
even  for  user  velocities  of  the  order  of  3,000  miles/hour,  due  to  the  low  transmis¬ 
sion  delays  induced  by  the  SALCS.  Thus,  for  the  same  end-to-end  propagation  delay, 
a,  our  simulations  are  identical  for  both  the  Ethernet  and  the  mobile  radio  LANs. 

Our  simulation  results  were  practically  identical,  when  the  SALCS  operational 
parameter  is  either  the  end-to-end  propagation  delay,  a,  or  the  mean  propagation 
delay,  a/2.  This  represents  a  positive  property  of  the  SALCS;  its  performance  is 
robust  regarding  the  choice  of  the  above  operational  parameter.  It  is  thus  suffi¬ 
cient  that  the  users  know  only  an  approximate  value  of  the  network  end-to-end  pro¬ 
pagation  delay,  for  the  SALCS  implementation.  In  figure  7,  we  plot  the  simulated 
expected  per  packet  delays,  D,  for  varying  Poisson  rates,  A,  together  with  the 
theoretically  computed  expected  per  packet  delays  induced  by  the  SALCS,  when  assump¬ 
tion  A1  holds  (table  2).  We  observe  that  the  SALCS  performance  is  better  in  the 
Ethernet  and  mobile  radio  environments,  than  in  a  star  network  environment  (where 
assumption  A1  holds).  This  performance  difference  increases,  as  the  network 
parameter,  y,  decreases,  reaching  its  maximum  for  y  *  0  (when  a  collision  of 
practically  zero  length  is  distinguished,  from  a  successful  transmission).  In 
addition,  in  the  Ethernet  and  mobile  radio  environments,  the  SALCS  out  performs  its 
synchronous  counterpart  ([18]),  as  concluded  by  comparing  figures  5  and  6  with 


figure  7. 
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7.  A  Refined  Algorithm 

As  we  pointed  out  in  section  3,  in  the  SALCS,  packets  do  not  use  the  distinction 
between  success  and  collision,  until  they  sense  the  first  idle  channel  state.  The 
question  that  arises  then  is:  If  such  a  distinction  is  possible  and  used,  will  it 
result  in  significant  performance  improvement?  To  answer  this  question,  we 
developed  and  analyzed  a  refined  algorithm,  named  GALCS,  adopting  again  assumptions 
Al  and  A2,  in  section  2,  The  GALCS  has  the  same  channel  sensing  characteristics  as 
the  SALCS,  and  it  is  again  implemented  independently  by  each  user  via  a  counter.  In 
addition  to  the  parameters,  a  and  y,  each  user  maintains  in  memory  seven  common 
parameters,  m^,  m2*  n' ,  pQ,  qQ,  r^,  and  however.  The  parameters,  m^,  n^,  and  n' 
are  positive  integers,  and  the  parameters,  pq,  qQ,  r^,  and  r^  are  probabilities; 
they  are  all  subject  to  optimization  for  throughput  maximization.  The  operation  of 
the  GALCS  is  described  by  parts.  A,  B,  and  C,  in  the  description  of  the  SALCS 
(section  3),  and  by  the  recursions  of  the  counter  value,  displayed  below. 


0,  then,  the  user  has  been  successfully  transmitted  at  t'. 
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If  at  (t')“  the  user  senses  the  channel  state  as  collision  and 

d)  rt  =  -1,  then, 

rt<=  0  5  with  probability  qQ. 

rti=  M  ;  1  £  M  £  m2,  with  probability  (l-qQ)  ni^^. 

e)  r  =  0,  and  the  user  acted  as  in  C.  b)  upon  sensing  the  collision,  then, 
r^i  =  0  ;  with  probability  rQ. 

r£i  =  M  ;  1  £  M  £  m2>  with  probability  r^-m^. 

rt'  **  m2+^  »  1  £  N  £  n'  ,  with  probability  (l-ro~r^)n'  * 

f)  rt,  =  0,  and  the  user  acted  as  in  C.  a)  upon  sensing  the  collision,  then. 


g)  rt  >_  1,  then 

rt,  «*  rt  +  m2  +  n' 

Rule  3:  After  time  t' ,  the  counter  value  is  updated  as  follows. 

a)  If  r  ,  >  1,  the  user  decrements  his  counter  value  by  one,  every  a  time 

units,  after  t',  to  a  value  not  less  than  0,  until  the  time  instant  t", 
that  he  first  senses  the  channel  busy  again.  Then,  he  maintains  the 
counter  value  rt„  unchanged,  until  he  first  senses  the  channel  idle  again 
He  then  repeats  rule  2. 

b)  If  r£  “  0  and  rt,  =1,  then, 

r  ■  -1  ;  for  t'  £  t  £  t*  +  a. 

rtl+a  =  0  ;  with  probability  rQ 

rtf+a  *  M  ;  1  £  M  £  m2>  with  probability  r^  nu*. 

r »  ,  *  m„  +  N  ;  1  <  N  <  n',  with  probability  (1-r  -r.)n'  1. 

b  *rCX  a  —  O  J. 

If  rt  *^>1,  the  user  decrements  his  counter  value  by  one,  every  a  time 
units,  after  t'+a,  to  a  value  not  less  than  0,  until  the  time  instant 
t",  that  he  first  senses  the  channel  busy  again.  Then  he  maintains  the 
value  rt««  unchanged,  until  he  first  senses  the  channel  idle  again.  He 
then  repeats  rule  2. 

We  observe  that  the  SALCS  evolves  from  the  GALCS,  via  the  following  substitutions 
in  the  latter:  m1*m2,*m-l,  pQ“q Q*m  \  r^r^O,  and  n'-n.  We  note  that  the  SALCS 
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represents  a  significant  operational  simplification  over  the  GALCS.  We  analyzed  the 
GALCS,  subject  to  assumption  Al.  We  found  that  the  improvement  in  performance  over 
the  SALCS  is  then  insignificant  (throughput  increase  in  the  third  decimal,  in  most 
cases),  at  the  expense  of  increased  operational  complexity.  The  above  statement  is 
also  true,  in  the  presence  of  the  mobile  radio  and  the  Ethernet-type  environments. 
Thus,  the  SALCS  is  practically  sufficient. 

8.  Conclusions.  Comparisons 

The  existing  classes  of  random  access  protocols  for  local  area  networks  are 
reflected  by  references  [3],  [5],  [12],  [18],  [22]  and  [23].  In  [3],  an  "Aloha- 
type"  protocol  is  studied.  For  the  Poisson  user  model  that  we  consider  in  this 
paper,  the  throughput  (the  maintainable  Poisson  intesity)  of  the  above  protocol  is 
zero.  The  protocol  proposed  in  [5]  belongs  to  the  "Aloha- type  with  retransmission 
control"  class  and  it  is  stable,  but  it  requires  full  channel  sensing  as  well  as 
knowledge  of  the  system  backlog,  and  it  is  synchronous  (slotted  channel) .  The 
protocol  in  [12]  is  basically  Gallager's  algorithm  [9],  with  different  cost  assign¬ 
ments  for  empty,  success,  and  collision  channel  states.  This  protocol  also  requires 
full  channel  sensing,  and  it  is  synchronous.  The  protocol  in  [22]  is  basically 
synchronous  and  it  requires  knowledge  of  the  system  backlog  (thus  complete  knowledge 
of  the  total  channel  history) . 

Considering  the  fact  that  in  most  local  area  networks  the  users  can  not  know 
the  total  channel  history  from  the  beginning  of  time  (before  they  enter  the  system) , 
the  only  sensible  approach  is  the  consideration  of  the  "limited  channel  sensing" 
class  of  protocols.  Those  protocols  only  require  that  each  user  sense  the  channel 
continuously,  from  the  time  he  generates  a  new  packet,  to  the  time  this  packet  is 
successfully  transmitted.  The  protocol  in  [18]  belongs  to  this  class;  it  requires, 
however,  slotted  channel  and  synchronous  transmission.  The  protocol  in  [23]  (ALCS) 
and  the  SALCS  in  this  paper  also  belong  to  the  limited  channel  sensing  class,  and 
they  address  the  case  where  users  can  not  distinguish  slot  boundaries,  due  to  lack 
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of  global  time  base.  The  ALCS  and  the  SALCS  are  thus  asynchronous,  they  are 
implemented  by  each  new  packet  independently,  and  they  induce  minimal  operational 
complexity.  The  difference  between  the  ALCS  and  the  SALCS  protocols  lies  in  the 
performance  superiority  of  the  latter.  In  star  networks,  where  the  propagation  delay 
between  any  two  users  remains  constant  and  unchanged,  the  SALCS  performs  worse  than 
its  synchronous  counterpart.  This  difference  in  performance  decreases  monotonically 
to  nil,  as  the  propagation  delay  decreases  to  the  order  of  thousandths  of  a  packet 
length,  and  it  is  then  compensated  by  the  elimination  of  the  synchronicity  requirement 
in  the  SALCS.  In  mobile  radio  and  Ethernet-type  environments,  the  SALCS  performs 
better  than  its  synchronous  counterpart,  since  it  then  models  the  real  environment 
closer  (see  also  [24]).  The  on  /  implementable  in  the  above  environments  existing 
stable  protocols  are  those  in  [18],  in  [23],  and  the  SALCS.  Among  them,  the  SALCS  is 
significantly  superior. 
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APPENDIX 

Proof  of  Proposition  1 

We  will  prove  equation  (7),  for  k  >  1.  The  proof  is  similar,  for  k  =  1.  Assume 

that  a  session  begins  at  time  0,  and  e ids  at  time  t'.  (Figure  8).  Let  us  also 

assume  that  the  first  subsession  in  the  session  is  of  multiplicity  k  >  1,  and  begins 

at  time  t  (case  b  Figure  8)  .  Let  us  denote  by  1^;  0  i  m+n-1,  the  number  of 

packets  contained  in  cell  //i  of  the  stack  at  time  t+&fot.  We  denote  by  t  the  time 

t+0+a,  and  by  t  the  time  t'  (for  uniformity  of  notation,  which  will  become  clear 
m+n 

in  the  sequel).  We  denote  by  t^  ;  1  _<  i  £  m+n-1  the  first  time,  during  the  session, 

when  the  1^  packets  occupy  cell  //I  of  the  stack,  while  cell  00  is  empty.  We  can  now 

define  the  random  variable  t.  .  ,  ;  1  <  i  <  ra+n  as  the  difference  between  the  r.v. 's 

i,i-l  -  - 

tJ  and  t.  ..  That  is,  t.  .  .  »  t.  -  t.  .  ;  1  <  i  <  m+n.  Furthermore,  L'(I  )  denotes 
i  i-1  i,i-l  1  i-l  —  —  o 

the  expected  value  of  t^  Q,  L"(Ii_1)  denotes  the  expected  value  of  ^  for 
2  _<  i  £  m,  and  L"'(I^  denotes  the  expected  value  of  t^  ^  ^  for  m  +  1  £  i  _<  m+n. 
Then  (see  also  figure  8)  we  have, 

m-1  m+n-1 

Lk  -  8+a  +  L'(Iq)  +  U(m-l)  •  +  X)  L"'^)  (A.l) 

i=l  i=m 

!1  for  m  >  0 
0  otherwise 

We  denote  by  H  .  0  <  i  <  m+n-1  the  number  of  new  arrivals  in  the  interval 

A,  i 

(t^,  tj+a) .  We  have  also  defined  as  L^  the  expected  length  of  a  subsession  of 
multiplicity  k  (k  _>  1) . 

We  now  have, 
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Consequently  (due  to  (A. 2),  (A. 3),  (A. A)) 
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Furthermore, 
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In  a  similar  way  we  can  show  that. 
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As  a  result. 
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Following  a  similar  procedure  we  can  show  that 
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Let  us  now  define  p(x,k)  =  e  X  yy  and  b^,(k,n)  =  (1  -  ~) 
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Then  (A.l),  (A. 5) ,  (A. 12)  and  (A. 13)  imply  that 
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When  k  =  1  we  take  (working  in  a  similar  way  as  when  k  >  1)  that 
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Proof  of  Theorem  2 


Part  (i).  For  convenience,  we  define  for  X  = 


°k«>  *  E  v  +  fk  •• k  i 1 


The  system  in  (7)  can  be  then  expressed  as  =  (^(X) .  Given  constants  b,  c  we 
construct  a  sequence,  Z^^  =  {z  such  that  z£^=bk-c  ;  k  _>  1.  After  straight 

forward  but  tedious  algebra  we  obtain, 
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o,  (Zvw)  =  x'w  +  g.  ;  k  >  1 


;  where. 
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kc  n  I 

We  now  define  the  sequence,  Z^n^  =  n  _>  1,  such  that,  z  =  0^(Z 

Since  O^O)  is  nonnegative  it  follows  that  for  every  k  _>  1,  the  sequence 


k  n>C 


nonnegative  and  nonincreasing,  if  b,  c  are  such  that. 


bk-c  >0  .V  k  >  1  and  e,  <  0  .  V  k  >  1 


Thus,  given  the  conditions  in  (A. 22),  the  following  limit  exists. 


lim  z,(n)  =  z,  <  bk-c  ;  k  >  1 
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The  sequence,  Z  =  {z  solves  system  (7).  Careful  examination  of  conditions  (A. 22) 

reveals  that  if  we  choose  b,  c,  such  that,  g^  =  =  0,  then  the  rest  of  conditions 

(A. 22)  are  satisfied  provided  that, 


f,  f1K 

r(X,a,8;m,n)  =  — | - f2b  >  0 

lc 


(A. 24) 


As  a  result,  X^Ca, 8;m,n)  is  the  unique  solution  of  the  equation  r(X,a,8;m,n)  =  0. 


Furthermore, 


,  r  r  f2c  fl  t  r  f2c  f lb  r  1  1 

b  =  l  f - - -  j  .  l  — - - f,L  J 


and. 
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c  =  —  +  f  b 
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(A. 25) 


(A. 26) 


The  construction  of  the  lower  bounds  in  theorem  2  is  similar  to  that  of  the  upper 
bounds,  and,  therefore,  it  is  omitted.  We  only  give  below,  the  expressions  for  b' 
and  c ' . 


-1 


lc 


oob 


(A. 27) 


where,  f  =  -dim  f,  ,  f  ,  =  £ira  f,  ,  ,  fm  =  £im  f, 

°°  k  °°b  , _  kb  “c  kc 


k-x» 


k-*» 


k-*°° 


and. 


(A. 28) 


Part  (ii)  In  this  part  we  show  that  the  expected  length,  L^,  induced  by  the  algorithm, 
coincides  with  the  solution  Z  =  (z  .  of  part  (i).  First  we  prove  the  following  lemma 


Lemma  A  The  solution  of  system,  -  0^(X),  is  unique  in  the  cJass  of  nonnegative 
sequences,  Y  =  which  satisfy  the  condition. 


£im  —j  =  0 
k-*»  k 


(A. 29] 


Proof  Assume  that  the  system,  x^  =  O^CX)  has  two  nonnegative  solutions  satisfying 
condition  (A. 29),  and  let  the  sequence  D  =  denote  their  difference.  Then, 

D  satisfies  (A. 29)  and  solves  the  following  system. 


dk  =  Yj  \,ydii  s  k  >  1 


For  u1  >  0  and  u^  k  +  u2k  +  u^  >  0  V  k  _>  1,  let  us  define  a  sequence  W  =  {w^J^ 
2 

such  that  w^  =  d^/Cu^k  +  u2k  +  u^)  ;  k  1. 

Then  the  sequence  W  solves  the  following  system. 


”k  '  E  \,y 
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(A. 30) 


_  2 

;  where  A^  ^  =  A^  ^/(u^k  +  u^k  +  u^)  ;  k  _>  1,  y  _>  1.  If  we  choose  u^ 
then  it  can  be  shown  that  there  exists  u^  >  0,  such  that. 


b,  u.  =  -c, 


A^  ^  <  1  ,  for  every  k  1 


(A. 31) 


Furthermore  we  know  that  w^  +  0  as  k  +  ®.  We  will  now  show  that  the  homogeneous 
system  (A. 30)  cannot  have  a  solution  tending  to  zero  and  different  from  zero.  Let 
us  denote  by  Q  >  0  the  exact  upper  bound  of  |wjJ  •  Since  w^  -*■  0,  we  will  have, 

beginning  with  some  k(k  _>  k^) ,  |w^|  <  Q/2,  and  therefore  it  must  be  possible  to  find 

a  k  <  k.  such  that  |w.  |  *  Q.  From  (A. 30),  for  k  =  k  ,  we  have. 
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Q  =  K  I  -  £  l\  I  l«l  1  Q  ]C  '"k  ,y 1 

O  .  o,y  y  , 

y=l  y=l 


(A. 32) 


From  (A. 31)  and  (A. 32),  we  have  that  Q  =  0,  and  therefore  w^  =  0,  for  every  k  _>  1. 
Thus  d^  =  0,  for  every  k  1,  and  the  proof  of  the  lemma  is  complete. 


The  solution  Z  of  part  (i)  clearly  satisfies  condition  (A. 29).  Thus,  from 
lemma  A,  we  have  that  Z  is  the  unique  nonnegative  solution  to  system  (7),  satisfying 
condition  (A. 29).  Next  we  use  arguments  similar  to  those  used  in  Theorem  2  of  [16] 
to  show  that  =  z^,  for  every  k  _>  1.  Let  us  consider  the  random  variables 
£k(x)  =  min()l^,T),  k  1  where  Jl^  is  the  subsession  length  of  multiplicity  k,  and  T 
is  a  real  number,  t  >  0.  Also  let  Lk(x)  =  E{£,^(t)};  k  ^  1.  It  can  be  easily  seen 
that,  for  every  k  _>  1: 


L^X^  <  Lk(T2)  if  0  <  xi  £  x2 


lim  L  (x)  =  L 


From  (A.  33)  and  from  the  dynamics  of  the  algorithm,  we  clearly  have  that, 

00 

Lk(x)  1  fk  +  ^  Ak, y  Ly(T^  *  for  every  k  >  1  or 
y=i 


Lk(T>  ■  fk  -  fk  +  E 


A.  L  (t)  ;  f.  _>  0  for  every  k  ^  1. 
k,y  y  k 


(A. 33) 


(A. 34) 


(A. 35) 


In  addition,  the  sequence  (Lk(t) )k>^  satisfies  condition  (A. 29)  for  every  T,  since 
L^Ct)  <_  X.  In  view  of  (A. 35),  this  fact  implies  that  Lk(x)  is  equal  to  the  difference 
of  the  nonnegative  solutions  zk  and  zk  of  system  (7)  that  respectively  correspond  to 
the  forcing  terms  f^  and  f^.  That  is,  Lk(x)  =  zk  ~  zk>  an<*  therefore  Lk(x)  £  zk  for 
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